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A method for the extraction of Lunar data and/or planetary 
features is provided. The feature extraction method can 
include one or more image processing techniques, including, 
but not limited to, a watershed segmentation and/or the gen- 
eralized Hough Transform. According to some embodiments, 
the feature extraction method can include extracting features, 
such as, small rocks. According to some embodiments, small 
rocks can be extracted by applying a watershed segmentation 
algorithm to the Canny gradient. According to some embodi- 
ments, applying a watershed segmentation algorithm to the 
Canny gradient can allow regions that appear as close con- 
tours in the gradient to be segmented. 
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AUTOMATIC EXTRACTION OF PLANETARY for crater detection can be divided into two main categories: 

IMAGE FEATURES supervised and unsupervised. The supervised methods 

require the input of an expert and generally use machine 
CROSS-REFERENCE TO RELATED learning concepts to train the algorithm to feature extraction. 

APPLICATION 5 Unsupervised methods are completely automatic and are gen- 

erally based on pattern recognition techniques. Different 
The present application claims a benefit from prior U.S. approaches have been proposed, based on template matching, 

Patent Application No. 61/179,774, filed May 20, 2009, texture analysis, neural networks, or a combination of these 

which is incorporated herein in its entirety by reference. techniques. Template matching has been described in A. 

to Flores-Mendez, “Crater marking and classification using 
STATEMENT OF GOVERNMENT INTEREST computer vision,” in Progress in Pattern Recognition, Speech 

and Image Analysis, vol. 2905, Lecture Notes in Computer 
The invention described herein was at least in-part made by Science. New York: Springer- Verlag, 2003 , pp . 7 9- 8 6, which 

an employee of the United States Government and may be is incorporated herein in its entirety by reference. Texture 

manufactured and used by or for the Government of the 15 analysis has been described in J. R. Kim, J.-P. Muller, S. van 
United States of America for governmental purposes without Gasselt, J. G. Morley, and G. Neukum, “Automated crater 

the payment of any royalties thereon or therefor. detection, a new tool for Mars cartography and chronology,” 

Photogramm. Eng. Remote Sensing, vol. 71, no. 10, pp. 

FIELD 13-22, 2000, which is incorporated herein in its entirety by 

20 reference. Neural networks have been described in A. A. 
The present invention relates generally to the field of auto- Smirnov, “Exploratory study of automated crater detection 

matically extracting features for the purpose of, for example, algorithm,” Boulder, Colo., 2002. Tech. Rep., which is incor- 

evaluating and comparing images, and particularly, lunar and / porated herein in its entirety by reference, 

or planetary images. Compared to Earth Science remote sensing data, lunar 

25 images usually present very low contrast and uneven illumi- 
BACKGROUND nation. The boundary of lunar features is not well defined, and 

it is therefore somewhat difficult to segment and characterize 
With the launch of several Lunar missions, such as Lunar lunar images. Also, because of uneven illumination, edges 
Reconnaissance Orbiter (LRO) and Chandrayaan-1 , a large extracted from lunar images do not form closed contours, and 
amount of Lunar images will be and are being acquired, and 30 post-processing needs to be done to link these edges . Further, 
will need to be analyzed. When registering or analyzing lunar because regions are difficult to characterize due to lack of 
data, significant features need to be extracted from the image contrast, if a method such as region growing is used, one level 
data. Planetary features, such as rocks, boulders, craters or of iteration is not sufficient to describe all the features. With 
ridges, are then used for applications such as: registering the large number ofnew lunar data that will be collected in the 
multi-temporal, multi-sensor, multi-view images; creating an 35 next few years, it is important to design an automated method 
obstacle distribution map for site selection or path planning to extract these features, and to perform tasks such as image 
purposes; or performing terrain categorization. Although registration. As such, an automated and robust feature extrac- 
many automatic feature extraction methods have been pro- tion method for lunar images is needed, 
posed and utilized for Earth remote sensing images, these 

methods are not always applicable to Lunar data that often 40 SUMMARY 

present low contrast and uneven illumination characteristics. 

The LRO is a NASA mission, aimed at creating a compre- According to various embodiments of the present teach- 

hensive atlas of the Moon features and resources to aid in the ings, a feature extraction method is provided, for example, for 
design of a lunar outpost and to prepare exploration and the extraction of Lunar data and/or planetary features. This 
scientific missions to the Moon. LRO is scheduled to spend at 45 feature extraction method can be used for extraction of fea- 
least one year in orbit collecting detailed information about tures, even when planetary images are blurry, noisy, present a 
the Moon and its environment. Different types of data will be lack of contrast and/or uneven illumination, and/or when the 
collected by LRO (and other Moon missions) at different represented images are not well-defined. According to some 

times, by different sensors, and from different view-points. embodiments, the feature extraction method can include 

Therefore, registration will be used to jointly exploit, inte- 50 computing an image gradient. The image gradient can be 
grate, or compare these different data, and feature extraction computed using a Canny edge detector, 
is the first step to not only image registration, but also any According to various embodiments, the feature extraction 

further analysis of these data. method can include one or more image processing tech- 

The identification of the features that are present on the niques, including, but not limited to, a watershed segmenta- 
planetary surface by a human expert is a time-consuming 55 tion and/or the generalized Hough Transform. According to 
endeavor. Therefore, a trustworthy automatic procedure to some embodiments, the feature extraction method can 

detect the position, structure, and dimension of each feature is include extracting features, such as, small rocks. According to 

highly desirable. This is a difficult task because limited data some embodiments, small rocks can be extracted by applying 

are available, the quality of the images is generally low (i.e., a watershed segmentation algorithm to the Canny gradient, 

it depends on illumination and surface properties), and the 60 According to some embodiments, applying a watershed seg- 
features that are present in the images can be barely visible mentation algorithm to the Canny gradient can allow regions 
due to erosion and exhibit different structures and variable that appear as close contours in the gradient to be segmented, 
sizes. According to some embodiments, the feature extraction 

Among typical features in Lunar- and planet- surface imag- method can include extracting features, such as, large rocks 
ery, craters play a primary role. The crater detection problem 65 and/or craters of elliptical shape. According to some embodi- 
has been widely addressed and different approaches have ments, large rocks and/or craters of elliptical shape can be 
been proposed in the literature. The image-based approaches extracted by applying a generalized Hough accumulator to 
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detect ellipses in the gradient image, and then applying the 
watershed segmentation using these ellipses as seed points. 

According to some embodiments, the feature extraction 
method can include detecting ridges. According to some 
embodiments, ridges can be detected by applying a standard 
Hough accumulator to detect straight lines in the gradient 
image. 

According to various embodiments, the feature extraction 
method can have many applications, which can include image 
registration. The feature extraction method can be general- 
ized to other planetary images as well as to Lunar images. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present teachings will be described with reference to 
the accompanying drawings. The drawings are intended to 
illustrate, not limit, the present teachings. 

FIG. 1 shows a flowchart for registering input and refer- 
ence images. 

FIG. 2A shows a partition of an infrared image. 

FIG. 2B shows an image of Canny Gradient of the image 
shown in FIG. 2A. 

FIG. 2C shows an image of segmented Canny Gradient. 

FIG. 2D shows an image of Canny Gradient with close 
contours extracted. 

FIG. 2E shows the resulting image after the segmentation 
(initialized with the maxima of the Hough accumulator) is 
applied. 

FIG. 2F shows a binary image in which crater boundaries 
are mapped. 

FIG. 3A shows an input image. 

FIG. 3B shows a reference image. 

FIG. 3C shows superimposed unregistered image bound- 
aries. 

FIG. 3D shows registered images by using a checkerboard 
representation. 

DETAILED DESCRIPTION OF THE PRESENT 
INVENTION 

According to various embodiments, a feature extraction 
method is provided that includes computing an image gradi- 
ent of a Lunar image by using an edge detector, detecting 
ridges in the image gradient by applying a Hough Transform 
to detect straight lines in the image gradient, and registering 
the image. In some embodiments, the method can include 
applying a watershed segmentation algorithm to the image 
gradient to extract object features of the image. The applying 
a watershed segmentation algorithm can include segmenting 
regions that appear as close contours in the gradient. In some 
embodiments, the edge detector can include a Canny edge 
detector. The method can further include detecting object 
features in the image, for example, rock features, crater fea- 
tures, ridges, or the like. In some embodiments, the method 
can further include preprocessing the image to smooth noise 
in the image prior to computing the image gradient. The 
registering can include registering the detected ridges with 
ridges shown in a reference image. 

According to various embodiments, a feature extraction 
method is provided that includes computing an image gradi- 
ent of an image by using an edge detector, applying a water- 
shed segmentation algorithm to the image gradient to extract 
first objects such as small rocks, applying a generalized 
Hough Transform to detect ellipses in the image gradient, and 
applying a watershed segmentation algorithm to the image 
gradient using the detected ellipses as seed points, for 
example, to extract objects and recesses of elliptical shape, 
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such as elliptical rocks and craters. The method can further 
include detecting ridges in the image gradient by applying a 
Hough Transform to detect straight lines in the image gradi- 
ent, and distinguishing one or more objects in the image from 
5 a background of the image. The edge detector can include a 
Canny edge detector. The image can include a Lunar image. 

In some embodiments, the method can further include 
preprocessing the image to smooth noise in the image prior to 
computing the image gradient. The method can also include 
to registering one or more of the detected objects with one or 
more objects in a reference image. 

According to various embodiments, the registering can 
include extracting features from both the input image and 
reference image and computing a geometric transformation 
15 between the image and the reference image based on these 
features. In some embodiments, the registering can include: 
extracting features from the reference image; extracting fea- 
tures from the input image; matching features between the 
input image and the reference image; estimating a transfor- 
20 mation model; transforming the input image with respect to 
the reference image; and generating a registered image. The 
generated registered image can be displayed, output, printed, 
written to a memory device, or a combination thereof. A 
processor can be used to carry out one or more of the extract- 
25 ing, the matching, the estimating, the transforming, and the 
generating. In some embodiments, a processor can be used to 
carry out two or more, three or more, or each of the extracting, 
the matching, the estimating, the transforming, and the gen- 
erating. In some embodiments, the method can further 
30 include printing out the registered image, displaying the reg- 
istered image, or the like. 

According to yet other embodiments of the present teach- 
ings, an image feature extraction system is provided that 
includes an edge detector, a processor, or both, for example, a 
35 processor having an edge detector software program installed 
therein. The processor can be configured to compute an image 
gradient of an image by using the edge detector. The proces- 
sor can be configured to apply a watershed segmentation 
algorithm to the image gradient to extract first objects. The 
40 processor can be configured to apply a generalized Hough 
Transform to detect ellipses in the image gradient. The pro- 
cessor can be configured to apply a watershed segmentation 
algorithm to the image gradient using the detected ellipses as 
seed points, to extract objects and recesses of elliptical shape. 
45 The processor can be configured to detect ridges in the image 
gradient by applying a Hough Transform to detect straight 
lines in the image gradient. In some embodiments, the pro- 
cessor can be configured to distinguish one or more objects in 
the image from a background of the image. 

50 In some embodiments, the system can include a Hough 
accumulator, a generalized Hough accumulator, or a combi- 
nation thereof. The system can include a Canny edge detector 
and the processor can be configured to compute the image 
gradient of the image using the Canny edge detector. The 
55 system can further include a camera configured to acquire the 
image, for example, configured to acquire a Lunar image. In 
some embodiments, the camera can be on board a spacecraft 
orbiting a planet, and the images can be of the surface of the 
planet. The processor can be configured to preprocess the 
60 image to smooth noise in the image. In some embodiments, 
the system can further include a preprocessor configured to 
preprocess the image to smooth noise in the image. 

According to various embodiments, a method for extrac- 
tion of Lunar features from Lunar images that show the sur- 
65 face of the moon and its structures, is provided. While Lunar 
features are specified, it is to be understood that the method 
can also be applied to other planetary features. According to 
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various embodiments, the method can enable automatic 
detection of structures that are present in a represented Lunar 
surface. The method can include extracting features and using 
the extracted features to register multi temporal, multisensor, 
and stereo -images. Different types of features can be present 5 
in the images. According to various embodiments, the fea- 
tures to be extracted can include rocks, craters, boulders, 
and/or ridges. Rocks are objects of small elliptical or circular 
shape, with no shadows. Craters have approximately ellipti- 
cal shape with shadows, due to their depth and uneven illu- to 
mination. Ridges appear like curves and straight lines in the 
images. The size, shape and position of the planetary images 
can be estimated by applying different methods, which are 
described herein. Once the features are extracted they can be 
applied to register image pairs, representing the same scene. 15 
As described in more detail herein, a genetic algorithm can be 
used to assess the transformation between the images to be 
registered. 

Given a pair of images to be registered (i.e., an input and a 
reference images), their features can be extracted and used to 20 
estimate the optimum transformation. Once the transforma- 
tion parameters are assessed, the transformation parameters 
can be used to register the input image with respect to the 
reference one. A flowchart of an exemplary technique is 
shown in FIG. 1 . As can be seen, features are extracted from 25 
both the input and the reference images (IN and REF, respec- 
tively). The extracted features are used in an optimization 
module, aimed at computing the optimum transformation. 
Once the transformation parameters are estimated, they can 
be used to transform IN with respect to REF and generate the 30 
registered image, REG. 

While extraction of features from lunar images is described 
herein, it should be understood that the methods described 
herein can apply to planetary features and planetary images as 
well. For example, the method for extraction of features can 35 
be used for planetary image registration, landing site selec- 
tion, hazard map creation and more generally for lunar and 
planetary terrain categorization. The method for extraction of 
features can further be used to supplement existing feature 
extraction methods already in use for Earth remote sensing, 40 
military and medical applications. In particular, the method 
can be useful for applications, such as, feature extraction in 
SAR images for military applications, and/or feature extrac- 
tion in medical images with low contrast (e.g., mammograms 
in MRIs). 45 

Feature Extraction 

According to various embodiments, the feature extraction 
method can include one or more image processing tech- 
niques. According to some embodiments, the image process- 
ing technique can include segmentation. The term “segmen- 50 
tation,” as used herein, can refer to a process of partitioning an 
image into multiple regions, for instance, in order to distin- 
guish objects from the background. According to some 
embodiments, segmentation can involve introducing a set of 
characteristic points that are related to the objects to be 55 
detected, automatically selected, and used as “seed points” to 
segment the images. Various approaches to segmentation can 
be used, for example, region growing, and/or, watershed 
transformation. 

According to various embodiments, the feature extraction 60 
method can include segmentation. According to some 
embodiments, the segmentation can include, for example, 
watershed segmentation, Hough Transform, and/or the Gen- 
eralized Hough Transform. 

According to various embodiments, watershed segmenta- 65 
tion can solve the problem of over- segmentation and can 
utilize good seed points from which regions can be extracted. 
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The Generalized Hough Transform can be used to extract 
ellipses to define these seed points. The watershed algorithm 
which can be used is described in S. Beucher, “The Watershed 
Transformation applied to Image Segmentation,” Scanning 
Microscopy International, 6, 1992, the contents of which are 
incorporated by reference in their entirety, herein. The water- 
shed algorithm is an automatic, robust, and fast method. The 
Generalized Hough Transform is described in S. Tsuji, F. 
Matsumoto, “Detection of Ellipses by a Modified Hough 
Transformation,” IEEE Trans, on Computers, 27(8), 1978, 
which is incorporated in its entirety herein by reference. 

According to various embodiments, before applying the 
feature extraction method, the images can be preprocessed. 
First, according to some embodiments, any noise that is 
present can be smoothed by applying a Gaussian filtering and 
a median filtering operation in cascade, as described, for 
example, in Shapiro, L. G. & Stockman, G. C: “Computer 
Vision”, Prentice Hall, 2001, which is incorporated herein in 
its entirety by reference. 

According to some embodiments, in order to detect the 
edges, the image gradient can be computed by using the 
Canny edge detector, described, for example, in J. Canny, “A 
Computational Approach to Edge Detection,” IEEE-Trans. 
Pattern Analysis & Machine Intelligence, 10(6), 1986, which 
is incorporated herein in its entirety by reference. According 
to some embodiments, both a gray scale image and a binary 
gradient can be obtained and used. According to some 
embodiments, as an intermediate result of this operation, an 
intensity gradient I can be generated, which can be a gray- 
scale image. According to some embodiments, by applying a 
non-maximum suppression algorithm followed by a hyster- 
esis thresholding to l g , a binary gradient image, \ b , showing 
the contours of the objects represented in the original image, 
can be obtained. 

According to various embodiments, in order to extract 
rocks, which can appear like close contours in the gradient 
image, the watershed algorithm can be applied to l b , in order 
to segment regions with close contours. According to some 
embodiments, all the area included within a close contour can 
be a “seed point-area”, and can be identified as a region. The 
result of this first step can be a binary image that shows the 
rock boundaries. 

While rocks can appear like close contours and can be 
easily detected, according to some embodiments, craters can 
have a more complex structure and, due to their depth and 
uneven illumination, can exhibit shadows. Their borders can 
be approximated with incomplete non-continuous elliptical 
curves. According to various embodiments, a generalized 
Hough accumulator, as described in S. Tsuji, F. Matsumoto, 
“Detection of Ellipses by a Modified Hough Transforma- 
tion,” IEEE Trans, on Computers, 27(8), 1 978, which is incor- 
porated in its entirety herein by reference, can be used to 
identify the seed points to detect these structures from I;,. For 
every pair of points that are detected as edge points in 1^ and 
exhibit opposite gradient direction, an accumulator, corre- 
sponding to the median point between them, can be incre- 
mented of a unit value. In other words, for each edge point 
pair with opposite gradient direction, an accumulator, corre- 
sponding to the median point between them, can be incre- 
mented. According to some embodiments, the maxima of the 
accumulator can be taken as centers of the ellipses. The 
parameters describing each ellipse centered in the detected 
maxima can then be assessed. According to some embodi- 
ments, a 3D accumulator can be used to estimate the two 
ellipse semi-axes and the rotation angle from all the pairs of 
points that contributed to the accumulator in the considered 
center. If the ellipse that has been generated truly corresponds 
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to a contour in the gradient image, its center can be used as a 
seed point for segmentation. According to some embodi- 
ments, starting from all the detected seed points, a watershed 
algorithm can be applied to I and the craters can be identified. 

As a result, a binary image that shows the crater boundaries 5 
can be obtained. According to some embodiments, a standard 
Hough accumulator can further be applied to detect straight 
lines in I . A standard Hough accumulator can be as 
described, for example, in R. Duda and R Hart, “Use of the 
Hough Transform to Detect Lines and Curves in Pictures,” 
Communications of the Association for Computing Machin- 
ery, 15, 1972, which is incorporated in its entirety herein by 
reference. 

According to various embodiments, ridges can be detected 1 5 
as unions of short linear segments. According to some 
embodiments, a binary image that represents the contours of 
all detected features, can be created and used in the registra- 
tion process. 

According to various embodiments, MATLAB can be used 20 
to develop computer code for implementing the methods 
described herein. In some embodiments, the watershed seg- 
mentation and the standard Hough Transform can use open 
source software. The generalized Hough transform can be 
loosely based on published algorithms. 25 

Registration 

According to various embodiments, once features are 
extracted, they can be applied to register image pairs repre- 
senting the same scene. Image registration, as used herein, 
can refer to a process of spatially aligning a pair of images. 30 
The pair of images can include an input image and a reference 
image. According to some embodiments, image registration 
can include a global optimization technique that is used to 
maximize the matching between features extracted from the 35 
images to be registered. Systems that evaluate and compare 
images, require image registration as an intermediate step. 
Images can be taken from different sensors, at different times, 
from different view -points. Registration steps can include 
feature detection, feature matching, transformation model 40 
estimation, and image transformation. Feature detection can 
include manual or automatic detection of salient or distinctive 
objects. Feature matching can include establishing corre- 
spondence between features. Transformation model estima- 
tion can include estimating a type of model and parameters of 45 
the mapping function. Image transformation can include 
transforming the input image. 

According to various embodiments, the feature extraction 
method can be applied to various registration approaches, 
including the basic approach, the point-based registration, 50 
and the region-based registration. The basic approach can 
include cross -correlation of gray -levels to obtain the measure 
of similarity between images. The point-based registration 
can include extracting and using peculiar points, for example, 
maxima or minima wavelet coefficients, to perform match- 55 
ing. 

According to various embodiments, the features described 
above can be extracted from the pair of images (input and 
reference) to be registered, and mapped into two binary 
images. Such binary images can be matched in order to com- 60 
pute the geometric transformation required to achieve the 
registration. According to some embodiments, the optimiza- 
tion scheme can include feeding binary images as inputs to an 
optimization module. The transformation matrix can be opti- 
mized, such that the evaluation can be performed using an 65 
objective function and the optimization can be achieved by 
applying a genetic algorithm. After the optimum matrix has 
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been estimated, it can be applied on the input image, for 
example, which can be translated and interpolated to obtain a 
final registered image. 

According to various embodiments, the feature extraction 
method can be applied to region-based image registration. 
According to some embodiments, the features described 
above can be extracted from the pair of images to be regis- 
tered, \ in and \ re p and mapped into two binary images, \ A and 
1^5, respectively. According to some embodiments, such 
binary images can then be matched to compute the geometric 
transformation required to achieve the registration. Accord- 
ing to some embodiments, registration can be based on a 
global optimization technique aimed at estimating the opti- 
mum parameters of an image transformation model. Accord- 
ing to some embodiments, the binary images, which represent 
the features of \ in and \ re p can be fed as inputs to an optimi- 
zation module. According to some embodiments, the trans- 
formation matrix can be optimized: its quality can be evalu- 
ated by an objective function and its optimization can be 
achieved by applying a genetic algorithm. According to some 
embodiments, after the optimum matrix has been estimated, it 
can be applied to one of the two images, which can be trans- 
lated and interpolated in order to obtain the final registered 
image. 

According to various embodiments, the problem can be 
formulated as determining a transformation T*such that, 
when T*is applied to the first image, l A , the best match with 
the second one, \ B is achieved. The match can be calculated as 
the correlation between the binary image extracted from the 
transformed input image and from the reference image. The 
objective function to be maximized is: 


M OM (T) = - V I A (T(x,y] I). ‘ 

n 

(x,y)\l g{x,y)± 0 

where MOM (measure of match) denotes the objective func- 
tion, T is the transformation for the x and y coordinates in the 
image plane, and n is the number of nonzero pixels of 1^. 
According to some embodiments, an affine transformation 
model, which exhibits six independent parameters, can be 
employed. The determination of the transformation param- 
eters can strongly depend on the objective function, as well as 
on the planetary images to be registered. According to some 
embodiments, where the function has multiple extremes, the 
most attractive search methods can be based on global opti- 
mization techniques. According to some embodiments, a 
genetic algorithm (GA) can be adopted, since it ensures, 
under mild assumptions, convergence to a global maximum 
of the adopted matching functional. The genetic algorithm 
can be as described in Z. Michalewicz, Genetic Algorithms+ 
Data Structures=Evolutional Programs, Springer Verlag, Ber- 
lin Heidelberg, third edition, 1999, which is incorporated 
herein in its entirety by reference. According to some embodi- 
ments, the six independent parameters of T can be defined 
over a wide range of values to achieve robustness. The aim of 
the GA can be to find the value for such parameters, which 
maximize the objective function. Given this, the final trans- 
formation matrix can be calculated by decoding the fittest 
individual of the last population, and the input image can be 
registered. 

For a description of other systems, methods, and/or appa- 
ratus, that can be implemented in accordance with various 
embodiments of the present teachings, reference is made to Y. 
Sawabe, T. Matsunaga, and S. Rokugawa, “Automated detec- 
tion and classification of lunar craters using multiple 
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approaches,” Adv. Space Res., vol. 37 , no. 1 , pp. 21 - 27 , 2006 ; 
J. Earl, A. Chicarro, C. Koeberl, P. G. Marchetti, and M. 
Milsen, ‘Automatic recognition of crater-like structures in 
terrestrial and planetary images,” in Proc. Lunar Planetary 
Sci. XXXVI , Houston, Tex., 2005 , Abs. No. 1319 ; and “Lunar 
reconnaissance Orbiter”, NASA Facts. Goddard Space Flight 
Center, National Aeronautics and Space Administration, 
http://lro.gsfc.nasa.gov/mission.btml; each of which is incor- 
porated herein in its entirety by reference. The present teach- 
ings are further illustrated with reference to the following 
examples which are intended to exemplify, not limit, the 
present teachings. 

EXAMPLES 

Experiments were carried out using data similar to that 
expected from LRO images but collected during the mission 
Mars Odyssey. Five bands visible images and ten bands ther- 
mal infrared images from the THEMIS (Thermal Emission 
Imaging System) instrument, with a resolution of 1 8 meters 
and 100 meters per pixel respectively, were used to test a 
method according to an embodiment of the present teachings. 

FIGS. 2A-2F shows results of a feature extraction method 
for a partition of the first band of an infrared image. The input 
image, shown in FIG. 2A, was first preprocessed, in order to 
smooth the noise. Canny gradient was applied to the 
smoothed image and results are shown in FIG. 2B. Subse- 
quently, a watershed algorithm was applied in order the 
extract the rocks. Segmentation results and the extracted rock 
boundaries are shown in FIGS. 2C and 2D, respectively. 
Finally, the generalized Hough transform was computed and 
a watershed segmentation was applied, which started the 
flooding process from the ellipse centers and allowed the 
craters to be detected. The segmentation results and the crater 
boundaries are shown in FIGS. 2E and 2F, respectively. 

To demonstrate the applicability of the method to registra- 
tion, two different non-registered bands of an infrared image 
were used. In order to show the results, the same partition of 
FIG. 2A was used. Results are shown in FIGS. 3A-3D. In 
FIGS. 3 A and 3B, the fourth and fifth bands were used. The 
rotation and translation between the two different bands are 
visible in FIG. 3C, in which two non-registered contour 
images are superimposed in a gray-scale representation. The 
contours extracted from the fourth band image are repre- 
sented in white whereas the fifth-band contours are shown in 
gray. Finally, the transformation parameters are estimated by 
the method and the co -registered images are shown in FIG. 
3D, by using a checkerboard filter, that is, each check of the 
board represents the registered input image and the reference 
image, alternately. The registration accuracy can be evaluated 
by looking at the continuity of the image features at the 
borders of the checks. The visual analysis of FIG. 3D shows 
that the registration performed well; craters and ridges appear 
continuous at the check borders, or points of overlap. 

Other embodiments of the present invention will be appar- 
ent to those skilled in the art from consideration of the present 
specification and practice of the present invention disclosed 
herein. It is intended that the present specification and 
examples be considered as exemplary only with the true 
scope and spirit of the invention being indicated by the fol- 
lowing claims and equivalents thereof. 

What is claimed is: 

1. A feature extraction method, the method comprising: 

computing an image gradient of a Lunar image by using an 
edge detector; 

detecting features in the image gradient by applying a 
Hough Transform to the image gradient; 
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at least one of mapping the features to generate a map, 
printing the features to generate a printout, and using the 
detected features to register the Lunar image and gener- 
ate a registered Lunar image; and 
5 applying a watershed segmentation algorithm to the image 

gradient and registering the Lunar image against an 
image of an object to extract an object feature of the 
Lunar image. 

2. The method of claim 1, wherein the features detected are 
ridges and the applying a Hough Transform comprises using 
a standard Hough accumulator. 

3. The method of claim 1, wherein the edge detector com- 
prises a Canny edge detector. 

4. The method of claim 1, wherein the method comprises 
using the detected features to register the image and generate 

15 a registered image. 

5. The method of claim 1, the method further comprising: 
preprocessing the Lunar image and the image of an object 

to smooth noise in the pair of images prior to computing 
the image gradient. 

20 6. The method of claim 1, wherein the registering com- 

prises registering the image to a reference image. 

7. A feature extraction method, the method comprising: 
computing an image gradient of an image by using an edge 

detector; 

25 applying a watershed segmentation algorithm to the image 

gradient to extract first objects; 
applying a generalized Hough Transform to detect ellipses 
in the image gradient; 

applying a watershed segmentation algorithm to the image 
30 gradient using the detected ellipses as seed points, to 
extract second objects and recesses of elliptical shape; 
detecting ridges in the image gradient by applying a stan- 
dard Hough Transform to detect straight lines in the 
image gradient; and 

35 distinguishing one or more second objects in the image 
from a background of the image. 

8. The method of claim 7, wherein the edge detector com- 
prises a Canny edge detector. 

9. The method of claim 7, wherein the image comprises a 
40 Lunar image. 

10. The method of claim 7, the method further comprising: 
preprocessing the image to smooth noise in the image prior 

to computing the image gradient. 

11. The method of claim 7, the method further comprising: 
45 registering the image with one or more reference images. 

12. The method of claim 11, wherein the registering com- 
prises extracting features from the reference image and com- 
puting a geometric transformation between the image and the 
reference image. 

50 13. The method of claim 11, wherein the registering com- 

prises: 

extracting features from the reference image; 
matching features between the image and the reference 
image; 

55 estimating a transformation model; 

transforming the image with respect to the reference 
image; and 

generating a registered image. 

14. The method of claim 13, wherein a processor is used to 
60 carry out one or more of the extracting, the matching, the 

estimating, the transforming, and the generating. 

15. The method of claim 13, wherein a processor is used to 
carry out each of the extracting, the matching, the estimating, 
the transforming, and the generating. 

65 16. The method of claim 13, the method further compris- 

ing: 

printing out the registered image. 
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17. The method of claim 13, the method further compris- 
ing: 

displaying the registered image. 

18. An image feature extraction system, the system com- 
prising: 

a processor configured to 

compute an image gradient of an image by using an edge 
detector, 

apply a watershed segmentation algorithm to the image 
gradient to extract first objects, 

apply a generalized Hough Transform to detect ellipses 
in the image gradient, 

apply a watershed segmentation algorithm to the image 
gradient using the detected ellipses as seed points, to 
extract second objects and recesses of elliptical shape, 
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detect ridges in the image gradient by applying a Hough 
Transform to detect straight lines in the image gradi- 
ent, and 

distinguish one or more second objects in the image 
5 from a background of the image. 

19. The system of claim 18, the system further comprising: 
a Canny edge detector, wherein the processor is configured 
to compute the image gradient of the image using the 
Canny edge detector. 

10 20. The system of claim 18, the system further comprising: 

a camera configured to acquire the image. 

21 . The system of claim 18, the system farther comprising: 
a preprocessor configured to preprocess the image to 
smooth noise in the image. 



